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A  siirple  instrument  is  described  that  meMures  the  intensity  of  solar 
radiation  in  four  spectral  regions.  Instruments  of  this  type  have  been  used 
to  observe  the  atmospheric  extinction  from  hi^  altitude  balloons  at  sunrise 
since  mid  196S.  Ihe  data  from  these  balloon  flints  are  described  and  the 
vertical  profiles  of  atmospheric  aerosols  deteimined  from  these  measurements 
are  presented,  IVidenoe  is  given  for  Hie  time  variation  of  aerosols  at  hi^ 
altitude. 


I .  Introduction 


In  the  past  decade  a  nunber  of  vertical  soundings  using  ballocn  borne  collecting 
devices  and  sampling  devices  tave  been  to  detemine  the  vertical  distribution, 
concentration  and  composition  of  the  atmospheric  aeroeols*  Also,  i'idirect  measurements 
have  been  made  of  the  aerosols  in  the  stratosphere  using  optical  techniques,  e.g. 
twili^t  sky  intensity  meeeurements,  searchli^t  and  laser  probing  and  solar 
extinction  measurements  from  rockets  and  satellites.  From  these  measurements 
it  has  long  been  realized  that  the  concentration  of  the  aerosols  varies  in  time. 

Because  of  the  oonplexity  of  the  collecting  and  sanq>ling  devices  and  the  difficulty 
of  reliably  observing  the  stratospheric  cierosols  from  the  surface,  it  has  been 
difficult  to  monitor  the  time  variations  of  this  important  atmospheric  constituent, 

A  small  simple  optical  instnsnent  vas  developed  that  measures  the  extinction 
of  the  solar  radiation  at  four  wavelengths  from  a  hi^  altitude  balloon.  From 
these  measurements  using  the  techniques  described  in  this  paper,  the  vertical 
concentration  times  the  effective  scattering  size  of  the  aerosols  can  easily  be 
determined.  Since  the  instrument  is  small  and  wei^  less  than  12  pounds  it  can 
routinely  be  carried  to  hi^  altitudes  on  100,000  ft^  balloons  enabling  the  study 
of  the  time  variations  of  the  aerosols, 

II,  Description  of  Equipment  and  Measurement, 

The  balloon  borne  instrvmant  consists  of  four  telescopes  that  look  in  the 
horizontal  direction  with  14  degree  x  14  degree  fields  each  centered  90  degrees 
in  azimuth  from  each  other.  This  asseirbly  is  rotated  beneath  the  ballocn  at 
a  constant  angular  speed  of  approxuretely  1  revolution  per  minute  as  the  ballocn 
floats  at  constant  altitude.  During  "Oris  time  the  rising  sun  traverses  the 
vertical  fields  of  the  telescopes.  The  intensity  of  tiie  solar  radiation  is 
measured  in  fotn?  spectral  regions  as  a  function  of  time  bv  the  telescopes  as 
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they  rotate.  The  intensities  are  meeisured  to  an  accuracy  of  1%. 

The  constructian  details  of  these  telescopes  are  illustrated  in  Figure  1« 
For  convenience  the  telescopes  will  be  referred  to  as  A)  B,  C  and  D  for  the 
remainder  of  this  report.  T^le  1  sunaariaas  liie  spectral  sensitivities  of 
the  telescopes. 


Thble  1 


T'elescope 

XC^A] 

AX  [OA] 

A 

3670 

300 

B 

4300 

350 

C 

5950 

400 

D 

9100 

1400 

(X  is  wavelength  of  maxiinLin  sensitivity  and  AX  is  the  full  width 
at  one  half  max  sensitivity. ) 

Telescopes  A,  B  and  C  eit5)loy  type  U29  vax3uum  phototubess  with  S-4  response 
as  detectors.  Telescope  D  uses  a  type  925  vacuum  phototube  vihich  hcis  a  S>1 
response.  For  each  telescope  the  vignetting  curve  is  measured  by  scanning 
over  the  svn  before  each  flij^t.  This  allows  one  to  make  a  oorrection  for 
any  deviation  from  a  flat  response  across  the  field.  This  correction  is 
*  less  than  4%. 

Figure  2  illustrates  the  exponentied  extinction  coefficient  per  atmosphere  (k) 
as  a  function  of  wavelength  produced  by  the  various  atmospheric  £d36orbers  under 
good  average  conditions  at  sea  level.  There  are  three  curves  illustrated  in 

4 

this  figure,  one  for  the  molecular  or  Raylei^  component  per  atmosphere  (k  «  1/X  ), 
one  for  the  aerosol  ocmpcnent  per  atmosphere  Ctiiis  curve  assumee  k  «  1/x  and 
is  normalized  to  C.  W,  Allen's  (1963)  value  at  lu)  and  the  third  c\irve  illustrates 
the  atmospheric  absorption  for  ozone  as  a  function  of  mvelength  assuming  3nn  of 
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ozcne  per  atmosphere.  Also  marieed  on  this  figure  eure  the  spectral  responses 
of  the  four  telescopes.  One  observes  that  telescope  A  is  priinarily  sensitive 
to  the  extinction  produced  by  the  molecular  or  Raylei^  8oattei?ing  of  the 
atmosphere.  Telescope  B  is  also  sensitive  to  the  ?aylei^  oon^xnent  and  to  some 
extent  the  extinction  produced  by  the  atmospheric  aierosols.  Telesocqoe  C  is 
centered  in  the  middle  of  the  Chappuis  bands  of  ozone.  Along  with  being 
sensitive  to  extinction  produced  by  the  molecular  and  aerosol  conqxnents  of 
the  atmosphere  it  is  sensitive  to  absorption  caused  by  the  Chappuis  bands. 
Telescope  D  is  primarily  sensitive  to  the  extinction  caused  in  the  atmosphere 
due  to  the  presence  of  the  atmospheric  aerosol  component. 

Figure  3  demonstrates  the  intensity  tiiat  one  would  expect  to  doserve 
as  a  function  of  atmospheric  path  for  each  of  the  four  telescopes  assuming  a 
pure  molecular  atmosphere. 

From  the  flight  of  April  12,  1967  Figure  4  hais  been  constructed  by  plotting 
the  observed  intensity  from  each  of  the  telescopes  as  a  function  of  the  calculated 
air  mass  (air  mass  calculations  will  be  described  in  the  next  secticn).  One 
ctoserves  the  similai?ity  of  the  cbta  from  telescopes  A  and  B  to  that  expected 
for  a  Raylei^  atmosphere  vhereas  the  data  of  telescopes  C  and  D  appear  different. 
The  difference  in  the  data  of  telescope  C  frcm  that  expected  for  a  Raylei^ 
atmosphere  is  a  reflection  of  the  extinction  produced  by  the  ozcne  ccnponent 
of  the  atmosphere.  The  inflection  in  the  data  from  telescope  D  betvgeen  5  and 
8  atmosbheres  is  produced  by  an  increased  ccnoentreticn  of  atmospheric 
aerosols  £dx}ve  the  tropopause  as  will  be  discussed  in  later  sections. 

The  data  from  this  series  of  measurements  from  telescope  C  will  be 
discussed  in  a  future  report  as  this  report  will  primarilv  conoern  itself  with 
the  atmospheric  aerosols. 
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III.  Calculation  of  Air 

order  to  plot  the  telesoope  data  as  a  fmction  of  air  inass«  <me  must 
calculate  the  atmospheric  air  mass  traversed  by  the  solar  ray  before  it  enters 
the  telescope. 

A  knotvledge  of  the  position  of  the  balloon  as  a  functicsn  of  time,  vhich 
is  determined  from  the  pressure  mecisured  using  an  OUand  cycle  on  the  balloon 
and  from  tracking  the  balloon  using  a  G.M.D. ,  enables  one  to  deteimne  by  the 
use  of  the  Air  Almanac  the  true  elevation  angle  (a)  of  the  sun  as  a  function 
of  time.  The  air  mass  can  then  be  calculated  knowing  the  elevation  angle  of 
the  sun  and  the  geometric  altitude  of  the  balloon  (hg).  The  geometry  of  the 
calculation  is  illustrated  in  Figure  S.  a  is  negative  for  elevation  angles 
below  the  horizontal  direction. 

Considering  first  positive  elevation  angles,  the  air  ness  traversed  by 
a  solar  ray  before  entering  the  telescope  at  altitvide  of  hg  when  the  sun  is 
at  elevation  angle  a  is  given  by: 

P(ho) 

A(a,hg)  =  A<a,0)  ,  for  a  5^  0  (1) 

where  P(h)  is  the  pressure  at  altitude  h  and  A(a,0}  is  the  nuirber  of  air 
masses  from  the  surface  at  elevation  angle  a  without  ihe  correction  for 
refraction.  Since  the  refraction  at  positive  elevation  angles  is  proportional 
to  the  density  at  the  balloon's  ed.titude,  ttie  correction  for  refraction  does 
not  significantly  influence  relation(l)  and  need  not  be  considered  for  positive 
elevation  angles.  The  surfaoe  air  mass  A(a,0)  has  been  measured  and  ced.culated 
for  various  model  atmospheres  by  a  nunter  of  authors.  (A  good  review  of  the 
determination  of  A(a,0)  has  been  presented  by  Roasnberg  (1966)).  For  the 
present  study  the  author  has  used  the  air  mss  tedsulated  by  Allen  (1963). 


For  negative  elevation  angles  one  can  calculate  the  air  mass  tieversed 
by  the  solar  beam  before  reaching  the  telescope  by  observing  that  at  the  tangent 
hei^t  of  the  ray,  hp,  the  air  mass  to  the  tangent  point  is  given  by: 

P(lV) 

A(0,h<p)  =  A(0,0)  t  (2) 


But,  the  toted  air  mass  at  elevation  angle  -a  is  just  twice  the  air  mass  to 
the  tangent  point  minus  the  air  mass  at  +a.  OR 


A{-a,hg) 


PClV) 

2A(0,0)  “  A(a,0) 


P(hR) 

TOJT"  * 


(3) 


In  relatians(2)  and  (3),hp  is  the  hei^t  of  the  tangent  ray  which  is 
affected  by  refraction.  Its  heij^t  can  be  calculated  by  an  iterative  procedure. 

From  Figure  S  it  my  be  seen  that  the  tangent  hei^t  of  a  ray  with  no 
refraction  is  given  by 

hp  =  hg  cos  a  -  R  (1-cos  a)  (4) 


where  R  is  the  radius  of  the  earth.  From  the  tangent  heij^t  for  a  straight  ray 
(no  refraction)  the  first  estirate  of  Ihe  pressure  of  -the  true  tangent  altitude 
is  determined,  and  thus  the  first  estimte  of  the  air  mass  (A’ )  to  the  tangent 
point  traversed  by  the  ray,  using  equation  (2).  The  refraction  of  the  tangent 
ray  to  the  tangent  point  is  then  calculated  as 

0  s  B  A'  - - -  (5) 

®  .965  +  .00  35  TChj.) 

vhere  T(hp)  is  the  temperature  at  altitude  hp  ^  degrees  centigrade.  In 
equation  (5), which  is  due  to  C.  W,  Allen  (1963),  0^  is  the  refraction  per  air 


TOSS  and  for  the  vavelengths  used  in  these  experiments  has  the  nvanerical  value 
1  min  of  arc  per  air  mas34 

Since  0  is  a  smll  angle  and  most  of  liie  refraction  occurs  near  the 
tangent  altitude  one  may  obtain  a  second  estimate  of  hj,  from  equation  (4)  bv 
substituting  for  a  the  apparent  elevation  angle  -a  ♦.0» 

Having  thus  calculated  the  tangent  hei^t,  corrected  for  refraction  to 
first  order,  relation  (2)  can  again  be  used  to  calculate  t^ie  air  mass  corrected 
for  refraction  to  first  order  and  then  the  calculation  of  the  tangent  heif^t 
can  be  repeated  and  the  iterative  procedure  can  be  continued. 

A  good  check  on  the  calculation  of  the  tangent  height  is  offered  by 
the  data  itself.  It  is  observed  that  in  flints  made  at  times  of  good  weather 
that  the  signal  of  telescope  D  is  first  observed  at  Idie  time  when  one  would 
just  predict  the  sun  had  come  over  the  limb  of  the  earth.  If  one  calculates 
the  tangent  height  at  the  time  tJie  first  signal  is  c^erved  using  the  above 
method  he  finds  it  is  only  a  fraction  of  a  km  indicating  ihe  oorrection 
for  refraction  has  been  applied  correctly  to  first  order. 

I\'.  Determination  of  Na(h)  from  Extinction  Data. 

The  nuRber  density  tines  the  effective  scattering  size  of  the  atmospheric 
aerosols  as  a  function  of  altitude,  Na(h),  oan  be  determined  frran  the  extinction 
data  from  telescope  D  by  considering  a  model  atmosphere  consisting  of  a  number 
of  layers  and  working  out  the  extinction  produced  in  each  of  "Ore  layers  as  the 
solar  rays  traverse  them. 

The  atmosphere  below  the  balloon  is  oonveniently  divided  into  a  nunber  of 
layers  by  the  tangient  heists  of  the  observed  rays  when  the  sun  is  below  the 
horizon.  It  is  oonvenient  to  index  the  layers  from  the  balloon  down  from  1  to 
k  and  assign  the  index  0  to  the  layer  above  the  balloon  of  thic)oiess  S  which  is 
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assumed  to  contain  the  aerosols  above  the  balloan*  The  tangent  hei^t,  hj^, 
can  be  calculated  by  the  method  outlined  in  the  previous  section  and  the  geometry 
is  illustrated  in  Figure  6.  The  path  lengths  of  the  ith  ray  in  the  kth  layer 
can  be  calculated  from  simple  geometry  and  it  is  found; 

=  [S  {S  ♦  2  (1^  ♦  R)  }]  (6) 

where  R  is  the  radius  of  the  earth, 

1/7 

Pjo  =  +  S  -  h^)  (h^  +  S  +  h.  +  2R)]  (7) 

1/2 

-  C(h  -  h.)  (h  +  h.  +  2R)] 
o  1  o  1 


and  for  k  >  0 

» 

1 

■■jk  =  2  I<\-l  -  N'  <\.l  *  ”  i‘  ^  '“ji  • 

i=k+l 

Invoking  the  Lambert-Beer  law  one  has  the  intensity  I  of  'Oie  li^t  ray 
after  traversing  a  path  length  P  when  the  incident  intensity  is  I'  as: 

I  =  I*  e-'^  P  (9) 

where  N  is  the  nunber  of  scattering  centers  per  cubic  volume  and  o  is  the 
effective  cross-section  of  the  average  scattering  center. 
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By  plotting  the  observed  intensity  from  telescope  D  corrected  for  the 

Rayleigh  caipcnent  as  a  function  of  air  mass  cne  can  extrapolate  to  find  the 

intensity  at  the  top  of  the  atmosphere,  Ig,  and  with  IJie  observed  intensity 

corrected  for  the  Rayleigh  ccmponent  in  ths  horiaontal  direction,  I  ,  one 

o 

can  inyert  relation  (9)  to  find 


by  using  relation  (6)  to  calculate  P^, 

Haying  detemined  (Ho)^,one  can  trace  the  first  ray  below  tee  horizontal 
through  the  atmosphere,  again  using  relation  (9)  edong  with  (7)  and  (8)  to  find 


lo^ 


■^10 


C)ne  can  continue  to  solye  for  tee  (Na)j^’s  using  the  meesured  intensities 
corrected  for  the  Rayleigh  corpcnent,  Ij^,  for  sucoessiye  lower  layers  and  in 
general  one  finds 


k-1 

'e  .  I 


loge  ®  (No).  P^. 

^  1  ki 

k  .1=0 


<Nah 


kk 


In  order  to  carry  out  the  above  cedculation  caie  must,  eis  has  been  pointed 
out,  assume  a  value  for  S.  In  the  present  work  S  was  first  assumed  to  be  the 
atmospheric  scale  height  and  then  after  the  calculation  was  completed  the 
assumption  was  made  that  o  vas  independent  of  eiltitude  and  the  scale  heif^t  of 
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the  aerosols  was  detennined.  This  scale  height  was  then  used  for  S  and  the 
oalculaticn  repeated.  Since  a  balloon  at  an  edtitude  of  32  km  is  above  almost 
all  of  the  aerosols  the  choice  of  S  does  not  critically  influence  the  results. 

As  a  test  of  the  consistency  of  this  method  in  one  of  tiie  balloon  flints 
of  this  series  telescope  C  was  replaced  with  a  telescope  T)  so  that  there  were 
tvo  identical  telescopes  looking  in  the  infrared  wavelength  on  the  same  fli^t. 
The  data  from  these  telescopes  were  analyzed  independentlv  and  are  plotted  in 
Figure  7.  Both  are  seen  to  give  the  same  result. 

V.  Data  from  Extinction  Flights. 

From  the  sunner  of  1965  to  early  1968  a  number  of  extinction  flints  were 
made  from  Minneapolis  and  one  fli^t  (774)  was  nade  from  the  Panama  Canal  Zone. 

The  data  from  these  flints  have  been  analyzed  using  the  method  outlined 
in  the  above  secticns  and  lkt(h)  Ms  been  detennined  at  the  time  of  each  of  the 
mBasurements .  Figures  8  through  29  illustrate  the  No(h)'s  determined  in  these 
measurements.  In  these  figures  the  (No)j^*s  have  been  plotted  vs.  the  tangent 
he.i^t  h.  for  the  k  th  layer.  The  curved  grid  lines  on  the  background  of  these 
figures  are  proportional  to  lines  of  constant  mixing  ratios  for  the  atmospheric 
aerosols  if  the  effective  scattering  cross-section  is  independent  of  altitude. 
The  numerical  values  on  the  lines ,  thus  interpreted,  are  proportional  to  the 
mixing  ratio. 

In  this  series  of  flights  the  hipest  concentratic»i  of  aerosols  above  the 
tropopause  at  Minneapolis  vas  observed  in  Fli^t  780  on  December  30,  1966,  The 
observed  No  in  this  layer  was  6  x  10"®  cm"^«  If  one  assumes  the  average  radius 
of  the  aerosols  responsible  for  the  extinction  was  .Sp,  one  calculates  the 
concentration  of  aerosols  in  this  layer  to  be  4/cm®. 
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In  Flight  774  from  Panana  a  hi^  Na(K>  6  x  10“®  cm"^)  was  c^Berved  in  the 
stratosphere  of  the  tropical  atmosphere.  This  flight  was  made  of  the  same  time 
that  J.  M.  Rosen  (1968)  made  direct  soundings  of  the  aerosols.  He  too  found 
a  hi^  concentration  of  aerosols  at  this  time  above  I^nama. 

VI.  Variability  of  the  Stratospheric  Aerosols. 

Study  of  these  figures  reveals  the  variability  of  the  eerosol  ccncentration 
in  the  stratosphere  and  at  an  altitude  iust  above  the  tropopause.  These 
variations  are  in  fact  systematic  in  time  as  is  evidenced  in  Figure  30.  In 
this  figure  the  extinction  coefficient  determined  fron  -Oie  (feta  of  the 
D  telescope  vdien  the  solar  ray  traversed  one  atmospheric  air  mass,  is  plotted 
above  the  data  of  each  fli^t.  As  has  been  pointed  out  in  earlier  sections, 

^IR*  ooef. ,  is  proportional  to  the  average  conoentraticn  of 

aerosols  along  the  path.  The  tangent  hei^t  of  a  path  throu^  the  atmosphere 
that  traverses  1  air  ness  is  about  26  km. 
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Constructicxi  details  of  telescope  used  in  extincticn  measurements. 
Exponential  extinction  coefficients  per  atmosphere  ais  a  function 
of  wavelength.  Spectral  sensitivities  of  the  telescopes  are  noted. 
Theoretical  intensity  as  a  function  of  atmospheric  path  for  the 
four  telescopes  in  a  pure  ^ylei^  atmosphere. 

Intensity  as  a  fxncticn  of  atmospheric  path  for  Flij^t  792. 

Geometry  of  the  air  mass  calculation. 

Geometry  in  a  layered  atmosphere  used  for  determining  No(h). 

Na<h)  for  the  fli^t  of  September  6,  1966.  Data  is  plotted  vhich 
was  determined  independently  from  two  different  telescopes. 

The  number  of  aerosols  times  the  effective  scattering  cross-section  No 
as  a  function  of  altitude. 

iDctincticn  coef.,  k_,  determined  for  the  26  km  layer  aa  a  function 
of  the  time  from  1965  to  1968. 
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